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analysis of the existing solar neutrino data. We obtained forbidden regions in the 



O-i. plane of the parameters Am^ and sin^ 29 in two cases: A) Without any restrictions 

r-| ! on the values of the solar neutrino fluxes from different reactions; B) With some 

restrictions that take into account the predictions of all the existing solar models. 
We show that the existing solar neutrino data allow to exclude rather large regions 
in the plane of the parameters Am^ and sin^ 29 (especially in case B) . In the 
second part of this report we present a general method for the analysis of solar 
neutrino data that can be applied to future solar neutrino experiments (SNO, 
Super-Kamiokande, Icarus) in which high energy ^B neutrinos will be detected. 
We show that these experiments will allow: 1) To reveal in a model independent 
way the presence of sterile neutrinos in the flux of solar neutrinos on the earth and 
to obtain lower bounds for the probability of transition of u^s into sterile states; 
2) To obtain directly from the experimental data the initial ^B z/e flux and the 
probability of z/g's to survive (if there are no transitions of z/g's into sterile states). 
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1 Introduction 

Solar neutrino experiments are very important for the investigation of neutrino mixing, as 
well as for the investigation of the sun. These experiments are sensitive to very small values 
of the neutrino mass difference squared Am^ = m^—ml (down to Am^ ^ 10^^'^ eV ) and to 
a wide region of mixing angles 6, including very small ^'s. On the other side, solar neutrino 
experiments allow us to detect neutrinos from different reactions of the thermonuclear solar 
cycles, including neutrinos from ^B decay, whose flux is about 10~^ of the total flux. The 
problem is that we cannot determine from the existing solar neutrino experiments [1-4] 
separately the values of neutrino masses and mixing angles and the values of the neutrino 
fluxes. Usually, to obtain information about the values of Am^ and sin^ 29, it is assumed 
that the values of the neutrino fluxes are given by the Standard Solar Model (SSM) [5-9]. 
It is well known, however, that the neutrino fluxes calculated in the framework of the SSM 
are subject to many sources of uncertainties, mainly due to a poor knowledge of some input 
parameters (especially nuclear cross sections and opacity). 
In this report we will present: 

1. The results of a solar model independent analysis of the existing solar neutrino 
data [10]. 

2. A model independent approach to future solar neutrino experiments in which solar 
neutrinos will be detected through the observation of CC, NC and neutrino-electron 
elastic scattering (ES) processes [11, 12]. 

We will show that future solar neutrino experiments (SNO [13], Super-Kamiokande [14], 
ICARUS [15]), in which high energy ^B neutrinos will be detected, will allow to answer 
in a model independent way the question whether there are transitions of solar z/g's into 
sterile states. If there are only active neutrinos Ug, z/^, Ur, in the flux of solar neutrinos 
on the earth, future solar neutrino experiments will allow to determine directly from the 
experimental data the initial flux of ^B z/g's and the probability of z/g's to survive as a 
function of neutrino energy. 

There exist at present data of three radiochemical solar neutrino experiments (Homes- 
take [1], GALLEX [2] and SAGE [3]) and the water Cherenkov direct counting Kamiokande 
experiment [4] . In the Homestake experiment solar neutrinos are detected through the ob- 
servation of the Pontecorvo-Davis reaction z/g + ^'''Cl — * e~ + ^^Ar, whose threshold is 
0.81 MeV. Thus pp neutrinos, which compose the main part of the solar neutrino flux, are 
not detected in this experiment and the main contribution to the Homestake event rate 
comes from ^B and ^Be (according to BP [5] 78% and 15%, respectively). In the GALLEX 
and SAGE experiments solar neutrinos are detected through the observation of the reac- 
tion z/e + ^^Ga — !► e~ -|- ^^Ge, whose threshold is 0.23 MeV. The main contributions to the 
GALLEX and SAGE event rates come from pp, ^Be and ^B (according to BP 54%, 27% 
and 8%, respectively). In the Kamiokande experiment solar neutrinos are detected through 
the observation of the process z/e — > z/e. The electron energy threshold in this experiment is 



Experiment Event Rate (SNU) SSM Predictions (SNU) 

BP TL CDF 



Homestake N^^y^ = 2.32 ± 0.23 8.0 ± 1.0 6.4 ± 1.4 7.8 

GALLEX N^l^ = 79 ± 10 ± 6 131.5+^ 123 ± 7 131 
SAGE Ar;f^L = 69±ll±6 

Kamiokande A^kam/^kam = 0-51 ± 0.04 ± 0.06 1 ± 0.14 0.8 ± 0.2 0.98 



Table 1: Data of solar neutrino experiments and rates predicted by BP [5], TL [6] and CDF [7]. 
The Kamiokande result is presented as the ratio of the observed rate to the rate predicted by BP. 

about 7 MeV. Thus only ^B neutrinos are detected in the Kamiokande experiment. In Ta- 
ble 1 the results of all four solar neutrino experiments are given. In the last three columns 
of Table 1 the values of the event rates predicted by Bahcall and Pinsonneault (BP) [5], 
Turck-Ghieze and Lopes (TL) [6] and Gastellani, Degl'Innocenti and Fiorentini (GDF) [7] 
are given. As it is seen from Table 1, the event rates measured in all four solar neutrino 
experiments are significantly less than the event rates predicted by the existing Standard 
Solar Models. 

Pontecorvo neutrino mixing (see, for example, Refs.[16, 17]) is apparently the most 
natural explanation of the possible "deficit" of solar neutrinos. It was shown in Refs. [18- 
24] that all existing solar neutrino data can be described by the resonant MSW mechanism 
[25] in the simplest case of mixing between two neutrino types^. The following two MSW 
solutions of the solar neutrino problem were found in the case of t'e^^'^(r) niixing: 

1. A small mixing angle solution: Am^ ~ 5 x lO^^eV^ and sin^ 2^ ~ 8 x 10~^ 

2. A large mixing angle solution: Am^ ~ 10^^ eV^ and sin^ 29 ~ 0.8. 

Thus, some indications in favor of neutrino mixing follow from the existing solar neu- 
trino data. We would like to emphasize, however, that these indications are model depen- 
dent: the analysis of the data is based on the assumption that the values of the neutrino 
fluxes from the different sources are given by the Standard Solar Models. 



^The data can be described also by vacuum oscillations [21] with a fine-tuning between Arn^ and the 
sun-earth distance. For the parameters Am^ and sin^ 29 the following values were obtained: Arn^ ~ 
8 X 10-" eV^ and sin^ 26 ~ 0.8. 



2 A model independent analysis of solar neutrino data 

In this section we present the results of a model independent analysis [10] of the data of the 
existing solar neutrino experiments. We will consider the simplest case of mixing between 
two types of active neutrinos [t^e^^nir)) and assume that the MSW mechanism takes place. 
Our analysis is based on the fact that the shapes of the spectra of neutrinos from the 
reactions of the thermonuclear cycle are known. These spectra are determined by the 
interactions responsible for the reactions and, as it was shown in Ref.[26], are negligibly 
affected by the conditions in the interior of the sun. Thus, the event rates of the solar 
neutrino experiments are determined by the values of the parameters Am"^ and sin^ 29 and 
by the values of the total neutrino fluxes (mainly pp, ^Be and ^B). We will consider the 
total neutrino fluxes as unknown parameters. From the existing solar neutrino data we 
cannot determine allowed regions of the values of the parameters Am^ and sin^ 29 and of 
the neutrino fluxes. Instead, we will obtain the regions of values of the parameters Am^ 
and sin^ 29 that are forbidden for all possible values of the initial neutrino fluxes (or for 
the initial neutrino fluxes constrained within wide limits chosen in such a way to include 
the predictions of all the existing Standard Solar Models). We will show that the existing 
solar neutrino data allow to exclude rather large regions of the values of the parameters 
Am"^ and sin^ 29. 

Let us write the initial spectrum of z/g from the source r [r = pp, pep, ^Be, *B, Hep, 
i^N, 1^0, i^F) in the form 

4>:^iE)=X^iE)<^^, (2.1) 

where E is the neutrino energy, X^{E) is a known function normalized by the condition 

/ X^(E) dE = 1 (see Ref. [5]) and $'' is the total initial neutrino flux from the source r. 

The integral event rate in any experiment a {a = HOM (Homestake), KAM (Kamiokande), 
GAL (GALLEX+SAGE)^) is given by the expression 

iV, = ^r;<|.'-. (2.2) 

r 

In the case of radiochemical experiments only solar z/g are detected. We have 

F; = / aa{E) X^{E) P..^..(E) dE , (2.3) 

with a = HOM, GAL. Here aa{E) is the cross section of the reaction z/g + ^^Cl -^ e^ +"^^Ar 
in the case of the Homestake experiment and the cross section of the reaction z/g + ^^Ga -^ 
e~ + ^^Ge, in the case of the GALLEX and SAGE experiments, Pi^^_>;y^(E) is the probability 
of Ue to survive and E^^^ is the threshold energy. For the calculation of the z/g survival 
probability we used the formula given in Ref. [27], which is valid for an exponentially 
decreasing electron density. 



^In our calculations we use the combined GALLEX-SAGE data: 74 ± 9 SNU. 



Source 


{E) 
(MeV) 


^LUM 

(MeV) 


$(BP) 

(cm^^sec^^) 


Smin 


Smax 


pp 


0.265 


13.10 


(6.00 ±0.004) X 10^° 


0.93 


1.07 


pep 


1.442 


11.92 


(1.43 ±0.02) X 10^ 


0.61 


1.29 


^Be 


0.813 


12.55 


(4.89 ± 0.29) X 10*^ 


0.46 


1.40 


SR 


6.710 


6.66 


(5.69 ±0.82) X 10^ 





1.43 


Hep 


9.625 


3.74 


1.23 X 10^ 


0.90 


1.13 


13N 


0.7067 


12.66 


(4.92 ± 0.84) X 10« 





1.51 


150 


0.9965 


12.37 


(4.26 ± 0.82) X 10^ 





1.58 


ITp 


0.9994 


12.37 


(5.39 ±0.86) X 10^ 





1.48 



Table 2: Solar neutrino fluxes (with la errors) predicted by BP; {E) is the average neutrino 
energy, Ilum = (3/2 — (i?), where Q = 26.73 MeV, and ^min and ^max determine the limits for 
the values of the total neutrino fluxes in case B. 

In the the Kamiokande experiment z/g as well as z/^ (and/or v^) are detected. We have 

YWm = I ^UE) X\E) P..^..(E) dE 

+ f a,^e(i^)X'(i5) $^P..^.,(E)dE, (2.4) 



=IM,T 



where a^^^eiE) is the cross section of the process z/^ e — ^ z/^ e (£ = e,/i), E^ is the recoil 
electron energy threshold and Pi^^^i^^(E) is the probability of the transition z/g ^ z/£ (£ = 
e,n). In our calculation we took into account the efficiency and the energy resolution of 
the Kamiokande detector [4]. The fluxes of neutrinos produced in the thermonuclear pp 
and CNO cycles must satisfy the following relation: 

iVLUM = 5^1i:uM<^^. (2.5) 

r 

HereA^LUM = Lq/AikP = (8.491±0.018) x 10^^ MeV cm'^ sec-\ where rf = 1 AU = 1.496 x 
10"^^ cm is the average sun-earth distance, Lq = (3.826 ± 0.008) x 10^^ erg sec ^""^ [28], is the 
luminosity of the sun and Vlum — Qf^ " {EY , where Q = 4 nip ± 2 mg — 'zn4He = 26.73 MeV 
and {Ey is the average energy of neutrinos from the source r. The values of {Ey and 
Vlum a-re given in Table 2. 

Our procedure for the analysis of the solar neutrino data is the following. At fixed 
values of the parameters Arri^ and sin^ 29 we calculate the x^ for all possible values of the 
neutrino fiuxes. For each value of the parameters Am^ and sin^ 29 and of the neutrino fiuxes 
we estimate the "goodness-of-fit" by calculating the confidence level (CL) corresponding 
to the calculated x^- Since we do not determine any parameter, the number of degrees 
of freedom of the x^ distribution is equal to the number of data points (i.e. four: three 
neutrino rates and the solar luminosity constraint). If all the confidence levels found for a 
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Figure 1 : Excluded regions in the sin^ 26- Am^ 
plane for MSW transitions due to t'e^f^rr) rnix- 
ing in case A. The region F is excluded at 95% 
CL within the solid line and at 95% CL within 
the dotted line. The allowed regions found with 
the BP neutrino fluxes are also shown (shaded 
areas) . 



Figure 2: Excluded regions in the sin^ 26- Am"^ 
plane for MSW transitions due to z/e-t'^(^) mix- 
ing in case B. The regions F are excluded at 
95% CL within the corresponding solid line and 
at 95% CL within the corresponding dotted line. 
The allowed regions found with the BP neutrino 
fluxes are also shown (shaded areas). 



given value of Am"^, sin^2^ and all possible values of the neutrino fluxes are smaller than 
a (we choose a = 0.1, 0.05, 0.01), then the corresponding point in the A?72^-sin^2^ plane 
is excluded at 100(1 — a)% CL. In this way we obtain the exclusion plots presented in 
Figs.l and 2. Let us notice that for the purpose of determination of the excluded regions 
in the parameter space our approach is the most conservative: any decrease of the number 
of degrees of freedom would increase the excluded regions. 

For the exclusion plot presented in Fig.l the only requirement was that all the total 
neutrino fluxes are positive. Let us call this case A. However, it is interesting and instructive 
to investigate how the forbidden regions in the Am'^-sm'^ 29 plane change if some limits on 
the allowed values of the neutrino fluxes are imposed. Thus we also considered the following 
case B: the different solar neutrino fluxes are constrained in the interval ^J^j^ $^(BP) < ^"^ < 
Cnax '^'^(BP), where $''(BP) are the BP values of the neutrino fluxes and the factors ^- 



r 
min 



and C,^s.x ^^^ chosen in such a way to include the predictions of the existing solar models 
[5-9]. The values of these factors are given in Table 2. We determined the minimum 
(maximum) values for the pp, pep, ^Be and Hep fluxes by subtracting (adding) 3 times the 
range of solar model predictions to the minimum (maximum) predicted flux (notices that 
this range is larger than the la error given by BP). Since it has been recently suggested 
[29] that the value of the astrophysical factor 5*17(0) could be significantly lower than that 
used in SSM calculations, we let the ^B flux to be arbitrarily small. Since the CNO fluxes 
have large uncertainties, we allow also them to be arbitrarily small. We determined the 
maximum values of the *B and CNO fluxes by adding 3 times the la error of BP to the BP 
average value. Let us emphasize that the limits on the allowed values of the neutrino fluxes 
which we imposed in case B are rather large. The excluded regions of the parameters Am^ 
and sin^ 29 in case B are presented in Fig. 2. As it can be seen from a comparison of this 
figure with Fig.l the excluded regions in case B are much larger than in the case where no 
limitation is imposed on the values of the neutrino fluxes. 

3 Are there sterile neutrinos in the flux of solar neu- 
trinos on the earth? 

The problem of existence of sterile neutrinos^ is very important for the theory beyond the 
Standard Model. Neutrino masses and mixing can be generated in the framework of the 
standard model if right-handed neutrino fields (singlets) together with left-handed doublets 
enter in the Yukawa interaction. In this case the total lepton charge L = L^ + L^ + L^ is 
conserved, neutrinos with definite masses are Dirac particles and transitions only between 
flavor neutrinos (z/^ -^ v^i with £, £' = e, /i, r) are allowed. In the models beyond the 
Standard Model transitions of active flavor neutrinos into sterile neutrino states become 
possible (see, for example, Refs.[16, 17]). So a discovery of such transitions will be a 
discovery of new physics. 

Here we will show that future solar neutrino experiments could allow to reveal the 
presence of sterile neutrinos in the solar neutrino flux on the earth independently on any 
assumption about the values of the initial neutrino fluxes [11, 12]. 

The solar neutrino experiments of the next generation (scheduled to start in 1996) are 
the SNO [13] and the Super-Kamiokande [14] experiments. In the SNO experiment solar 
neutrinos will be detected through the observation of three different processes: 

1. The CC process 

Ve + d ^ <^~ +P + P', (3.1) 

2. The NC process 

v + d^v + p + n; (3.2) 



•^Sterile neutrinos arc quanta of right-handed neutrino fields. In the Dirac-Majorana mixing scheme, 
these fields are mixed with the left-handed fields (see, for example, Refs.[16, 17]). Sterile neutrinos do not 
interact with matter via standard CC and NC interactions. 



3. The elastic scattering (ES) process 

z/ + e^ ^ z/ + e^ . (3.3) 

In the Super-Kamiokande (S-K) experiment solar neutrinos will be detected through 
the observation of the process (3.3) with an event rate about 30 times larger than in the 
current Kamiokande experiment. In both the SNO and S-K experiments, due to the high 
energy thresholds (~ 6 MeV for the CC process, 2.2 MeV for the NC process and ~ 5 MeV 
for the ES process), only neutrinos coming from *B decay will be detected. The energy 
spectrum of the initial *B z/g's is given by 

4>lSE)=X{E)^- (3.4) 

The function X{E) is the normalized neutrino spectrum from the decay ^B -^ ^Be + 
e"*" + z/e, which is determined by the phase space factor (small corrections due to forbidden 
transitions where calculated in Ref.[30]). The factor $ in Eq.(3.4) is the total flux of initial 
®B solar z/e's. 

Consider first the NC process (3.2). Using v^-v^-v^ universality of NC for the total 
NC event rate in the SNO experiment A^^*-^, we have 



iV^^= / o,^{E)X{E) V P,^^,,(E)dE<l> 



(3.5) 



where a^diE) is the cross section for the process u d ^ unp, E^ is the threshold neutrino 
energy and Puc^ue(E) is the probability of transition of solar z/e's into ug {i = e,fj,,T). It 
is useful to introduce the average total probability of transitions of z/g into other active 
neutrinos: 

L f aUE)X{E) V P.^^,,(E)dE, (3.6) 




Y 



where 



X,rf= / a,d{E)X{E)dE. (3.7) 



Ef^^ 



Using the results of a recent calculation of the cross-section a^diE) [31] we obtained X,jd ■ 
4.72 X lO-'^^cm^. From Eqs.(3.5) and (3.6) we get 

\£=e,ti,T I ^Q 

In the general case of transitions of z/g's into active as well as into sterile neutrinos 



£=e,fi,T 



P.,^.,(E) = 1-P.^^,,(E), (3.9) 



where P,y^^,yg(E) is the total probabihty of transition of z/e's into all possible sterile states^. 
From Eqs.(3.8) and (3.9) it follows that 

AtNC 

1 - (P.e^.s)Nc = x:;^ ■ (3-10) 

Thus, from the measurement of the NC event rate A^'^^'^ it is impossible to reach any 
conclusions about transitions of solar z/g's into sterile states without assumptions about 
the value of the total flux $. However, if solar ^B neutrinos are detected not only through 
the observation of NC but also through the observation of the ES and CC processes the 
problem of existence of sterile neutrinos could be solved in a completely model independent 
way. 

Let us consider the ES process (3.3). The total number of ES events is equal to 

iV^' = / crUE) P..^..(E) X(E) dE$+ / a,^e(E) J^ P-c-..(E) X(E) dE$ . (3.11) 

From Eq.(3.11) we have 

S^'= / ^uAE) E P.e^.,(E)X(E)dE<|.. (3.12) 

Jet .±^. 



Here 

J.ES ^ jyES 



/ {a,^,{E) - a^^^iE)) <P,XE) dE , (3.13) 

Jew 

where (f>ue{E) = P^^__»,^^(E) X(E) $ is the flux of z/g on the earth. The quantity S^^ can be 
obtained from the data of the SNO and S-K experiments. In fact, A^^^ will be measured 
in both experiments. In the SNO experiment the spectrum of the electrons in the CC 
process (3.1) will be measured and the spectrum of solar z/g on the earth, (f)i,^{E), will be 
determined. From Eq.(3.12) we obtain the relation 

which is similar in form to the relation (3.8). The quantities in the relation (3.14) are 
determined as follows: 

X,,e= / a,^,{E)X{E)dE (3.15) 

and 

E P'^e^'^O ^J^/es^^"^^^^^^^^ 5Z P.e^..(E)dE. (3.16) 



''Wc assume that neutrinos are stable particles. For a discussion of neutrino instability see Ref.[32] and 
references therein. 



For E^ = 5.94 MeV (which corresponds to a kinetic energy threshold Tth = 4.5 MeV for 
the electrons in the CC process) we have X^, g = 2.08 x 10~^^cm^. From Eqs.(3.8) and 
(3.14) we have 



where 

■■^Nc — ^^ atNC ■ yo.ia 






Let us stress that in the ratio R^|; only measurable and known quantities enter (the flux 
$ cancels in the ratio). From Eq.(3.17) it is evident that if 

Rnc ^ 1 (3.19) 

it would mean that there are transitions of solar z/g's into sterile states. 

In the case R^|; = 1 no conclusion about sterile neutrinos can be reached. In fact, 
the ratio R§|. is equal to one if {Pu^^ug)^^^ = {Pu^^ug)^^. This relation is satisfied at 
Puc->usO^) = as well as at Pu^-,i,g(E) = const ^ 0. Thus, if R^c 7^ 1 it would mean not 
only that sterile neutrinos exist, but also that the probability of the transition of solar z/g's 
into sterile states depends on neutrino energy. 

If the inequality (3.19) takes place, it is possible to obtain lower bounds for the average 
values of the probability of transition of z/g's into sterile states. In fact, we have 




E P--- ) ^ T^^ f^ = « • (3-20) 




From this inequality it follows that 

{Pv^~.us)es > 1 ~ Rnc • (3-21) 

Thus, if the inequality 

K% < 1 (3.22) 

is satisfied, from Eq.(3.21) we obtain a non-zero lower bound for the average probability 

{Puc^us/ES- 

In the case of transitions of solar z/g's into sterile states the initial ^B z/g flux cannot 
be determined from the experimental data. However, in this case we can obtain from the 
data a lower bound for this flux. In fact, from Eqs.(3.14), (3.18) and (3.20) we have 

^ES ArNC 






^ > -T^-t;Wc = ^^ ■ (3-23) 



MSW 4.5 X 10~6 7.0 X lO^^ 0.57 0.71 

VACUUM OSC. 6.3 x lO"" 0.85 0.56 0.50 



Table 3: Results of the calculation of [Pj^e^^^elmax ^^'^ -'^cc ^^ ^^^ model with fg-fs mixing. The 
values of Am? and sin^ 26 used are also given. These values were obtained from the analysis of 
the existing experimental data (Ref.[20] for the MSW transitions and Ref.[21] for the vacuum 
oscillations). 

Further, we have 




(3.24) 



From Eq.(3.24) it follows that if the inequality 

Rfc > 1 (3.25) 

is satisfied, then for the average probabihty (Puc^us)j^c ^^^ ^°^ ^^^ total flux $ we have 
the following lower bounds: 

atNC yES 

$ > . = —— . (3.27) 

^ud (%cj ''^''' 

It is instructive to calculate the measurable ratio R^^ (and other ratios that we will 
consider later) in some model. We considered the simplest model with z/g^^s mixing. The 
parameters of the model ^rn? and sin^ 2^, which are given in Table 3, were obtained from 
the fit of the existing solar neutrino data under the assumption of MSW transitions or 
vacuum oscillations (assuming the SSM neutrino fluxes). The ratio R^|. as a function of 
the recoil electron threshold energy is depicted in Fig. 3. From Fig.3 it can be seen that 
in the case of MSW transitions the ratio R^|; is bigger than one for all the considered 
threshold energies. In both the MSW and vacuum oscillation cases the ratio R^|; increases 
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Figure 3: Results of the calculation of the ratio R^c (^^^ Eq.(3.18)) in the model with Ve~v<s, 
mixing (Tth is the kinetic energy threshold in the ES process). The curves correspond to MSW 
transitions and vacuum oscillations. The values of the parameters Atti^ and sin 2Q used in the 
calculation are given in Table 3. 



with Tth- From Fig. 3 it follows that a high threshold energy is preferable for revealing the 
existence of sterile neutrinos. 

A detailed investigation of the spectrum of the recoil electrons in the process ve — > ve 
will be carried out in the S-K and SNO experiments. We will discuss now what additional 
model independent tests of the existence of sterile neutrinos can be performed when the 
recoil electron spectrum will be available. The spectrum of recoil electrons is given by 



n^«(T) 



£m{T) 



dT 



■(E,T)P.^^,^(E)X(E)dE$ 



+ 



dcr„ 



(i?,T) J]P,^^,,(E)X(E)dE$. 



^-(T) dT ^^^^^ 



(3.28) 



Here T is the kinetic energy of the recoil electrons, -Em(T) = T ( 1 + y 1 + 2 me/T j /2 and 

— -^(-E,T) is the differential cross section of the process v^e -^ v^e (i = e,/i). With the 
dT 



help of Eq.(3.28) we get the following relation: 



Uc^Ul 



=e,ii,T 



ES;T 



S^^(T) 
*X,^e(T) • 



(3.29) 
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Figure 4: Plot of the functions X,y^eC^) and Figure 5: Results of the calculation of the ra- 
X^^e(T) defined in Eqs.(3.30) and (4.11), re- tio R^S(T) (see Eq.(3.34)) in the model with 
spectively. The depicted range for the kinetic i^e~J^s rtiixing (T is the kinetic energy of the re- 
energy T of the recoil electrons in the ES process coil electrons in the ES process) . The curves 
will be explored by SNO with Tth = 4.5 MeV. correspond to MSW transitions and vacuum os- 
cillations. The values of the parameters Am^ 
and sin^ 26 used in the calculation are given in 
Table 3. 



Here 



x.M) 



Em(T) 



dT 



(E,T)X(E)dE 



(3.30) 



is a known function, which is plotted in Fig.4. Other quantities in Eq.(3.29) are determined 
as follows: 



E'=S(T) = n^S(T) 



Em(T) 



^^!^(E,T)-^^^(E,T) 
dT ^ ' ^ dT ^ ' ^ 



0,,(E)dE (3.3i; 



and 



Ep. 






do-. 



Xu^eO^) JEmm dT 



i(i?,T)X(E) J2 P.e-..,(E)dE. (3.32) 

£=e,fj,,T 



Let us stress that to determine the quantity S^^(T) it is necessary to know the recoil 
electron energy spectrum n^^(T) as well as the spectrum of z/g on the earth (f)„^{E). 

Let us consider now the relation (3.29). If the quantity S^^(T)/Xj, c(T) depends on the 
energy T, it would mean that transitions of solar z/g's into sterile neutrinos take place (if 
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there are no such transitions, then 



Y, P.e-.., ) = 1 and S^«(T)/X,^e(T) = const). 

In order to derive a lower bound for the average value of the probability of transition of 
z/g's into sterile states we will use the following inequality: 



Ep. 






p.. 






< 






Here 



R^^(T) 






R^^(T) 



(3.33) 



(3.34) 



[S"'(T)/X..e(T)]_ 

is a measurable quantity and the subscript max indicates the maximum value in the ex- 
plored energy range. From Eq.(3.33) it follows that 

(P,,^,3).g^T > 1 - R^S(T) . 



(3.35) 



For the total initial flux of ^B neutrinos, from Eqs.(3.29) and (3.33) we obtain the 
following inequality: 

$> — ^ . (3.36) 

The results of the calculation of the ratio R^^(T) in the model with fe^^s mixing is 
presented in Fig. 5. From this figure it is seen that a detailed investigation of the ratio 
R^^(T) near the threshold could be a promising way to search for z/g — > z/g transitions. 

We will obtain now other inequalities the test of which could allow to obtain model 
independent informations about the existence of sterile neutrinos. We have 



Ep. 



= e,tJ.,T I J.grp 



< 



Here a = NC, ES and the ratios 






Rf(T) 



Rn'c(T) = 



S^^(T)X 



ud 



X,^e(T)N^^' 

S^'(T)X,„ 



X.,c(T) S 



ES 



(3.37) 



(3.38) 
(3.39) 
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are measurable quantities. 
If the inequality 

Rnc(T) < 1 (3.40) 

takes place in some region of T, it would mean that there are transitions of solar z/g's 
into sterile states. From Eqs.(3.29) and (3.37) we find the following lower bounds for 
(P;y^^^3)gg.rp and for the initial total Ue flux: 

(P..^.s)es;T>1-Rnc(T), (3.41) 

AtNC 
$ > -^ . (3.42) 

Analogously, if in some region of (T) the inequality 

ReKT) < 1 (3.43) 

is satisfied, we have the following lower bounds for (Pv^^vs)-es-t ^^'^ *^- 

(P^c-<^s)ES;T>l-REi(T), (3.44) 

VES 

$ > -^ . (3.45) 



Further, we have 






\i=e,^i,T 



where a = NC, ES. Thus, sterile neutrinos exist if in some region of the variable T the 
inequality 

Rf (T) > 1 (3.47) 

is satisfied. For the averaged probability of the transition of Ue into us and for the total 
flux we have 

(P..-..3).>l-(Rf(T)):l, (3.48) 

~ \X^Mc(T)y„^^ 

In Figs. 6 and 7 we plotted the ratios R§|;(T) and Rfg(T) calculated in the model with 
z/g-z/s mixing (the parameters of the model are given in Table 3). It can be seen from these 
figures that in the model under consideration the ratios R§c(T) and Rfg(T) are larger 
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Figure 6: Results of the calculation of the ratio Figure 7: Results of the calculation of the ratio 
R^q(T) (see Eq.(3.38)) in the model with z^e^J^s R'EsC^) (^^^ Eq.(3.39)) in the model with J^e^J^s 
mixing. mixing. 



than one in a large part of the kinematical region. The deviation of R^|;(T) and RfKT) 
from one is larger in the case of vacuum oscillations than in the MSW case. 

We have derived several different inequalities whose test could allow to reveal the pres- 
ence of sterile neutrinos in the solar neutrino flux on the earth in a model independent way. 
Additional inequalities of this type can be obtained from the knowledge of the spectrum 
of solar neutrinos on the earth (\)yJ^E). The relation 



= (E) 



0.e(E) 



$X(E) 
and the relations (3.8), (3.14) and (3.29) have a similar structure. We have 



(3.50) 



Ep 



Vc^Vl 



< 



where a = NC, ES and the ratios 



R; 



NC 
CO 



Ep 



Uc^l^l 



=e,ii,T 



[P. 



•c^i'cJrnax 



N^^/X, 



ud 



pES 



. S^VX.^e 

15 



pa 
-■^CC ' 



(3.51) 



(3.52) 
(3.53) 
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Figure 8: Results of the calculation of the ra- Figure 9: Results of the calculation of the ratio 

tio RqI; (see Eq.(3.53)) in the model with fg- Rq^(T) (see Eq.(3.57)) in the model with Ue-vs 

fs mixing (Tth is the kinetic energy threshold mixing (T is the kinetic energy of the recoil elec- 

of the recoil electrons in the ES process). The trons in the ES process). The curves correspond 

curves correspond to MSW transitions and vac- to MSW transitions and vacuum oscillations, 
uum oscillations. 



are measurable quantities. The quantity [0i/e/-^]max ^^ Eq.(3.52) and (3.53) is the maximal 
value of the function (f),y^{E) / X (E) in the explored energy range. From Eq.(3.51), for the 
averaged values of the probability Pu^^us and for the total flux $ we have 






r; 



cc 



(a = NC, ES) , 



(3.54) 
(3.55) 



The results of our calculations of the ratio R^q in the model with z/g^t's niixing are 
given in Table 3, whereas the results of our calculations of the ratio R^^ as a function of 
the kinetic energy threshold of the recoil electrons in the ES process, Tth, are depicted in 
Fig. 8 (the kinetic energy threshold of the electrons in the CC process is assumed to be 
4.5 MeV). It can be seen from Fig. 8 that in the model under consideration the ratio R^c 
depends strongly on Tth (especially in the case of vacuum oscillations). It is preferable to 
search for z/g -^ ug transitions at relatively small thresholds. 

Finally, from Eqs.(3.29) and (3.50) it follows 



\£=e,fi,T 



Vc^Ul 



< Rcc(T) 



(3.56) 



ES;T 
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where the ratio 



Rc'c(T) - ^ ^tr'L (3.57) 



X,^e(T) [0.JX] 

is a measurable quantity. From Eqs.(3.29) and (3.56) we have 



(P..^.s)es;T>1-Rcc(T), (3.58) 

$ > ^ . (3.59) 

J max 

The resuhs of the calculations of the ratio Rcc(T) in the model with Ue^^s rnixing are 
presented in Fig. 9. In this model there are large deviations from one of the ratio Rcc(T) 
for values of T close to the threshold (especially in the case of vacuum oscillations). 

In conclusion, we would like to stress that all the inequalities discussed here are based 
on the conservation law (3.9). The concrete mechanism for the transition of z/g's into 
sterile states is not important for the tests proposed here. Thus, there could be not only 
'^e ~^ ^s transitions due to Dirac-Majorana neutrino mixing, but also u^ —^ vo, transitions 
due to neutrino Dirac magnetic moments [33]. In the latter case the measurable ratios 
R^c, R^^(T), . . .will depend periodically on time [12, 34]. 

4 The case of absence of transitions of solar v^^ into 
sterile states 

If the test of all the inequalities obtained in the previous section will not reveal the presence 
of sterile neutrinos in the solar neutrino flux on the earth, it will be natural to consider the 
possibility of a model independent treatment of solar neutrino data under the assumption 
that there are no transitions of solar z/p's into sterile states^. In this case 




Uc^Ul 



[a = NC, ES) , (4.1) 



u.^u, I = 1 , (4.2) 

and the SNO and S-K data will allow to determine the initial flux of ®B z/g's with three 
methods: 

1. By means of the measurement of the NC event rate A^'^^'^. In fact, from Eqs.(3.8) and 

(4.1) we have 

AtNC 

* = ^^ , (4.3) 

where X^^ is given by Eq.(3.7). 



''Let us stress that all the inequalities considered in the previous section will not reveal the presence of 
sterile neutrinos if P^^^^s(E) = const. 
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2. By means of the measurement of the total number of ES events A^^^ and the electron 
spectrum in the CC process (from which the spectrum of v^ on the earth, (f)i,^{E), 
will be determined). From Eqs.(3.14) and (4.1) we have 



VES 

$ = -^ , (4.4) 

where S^^ is given by Eq.(3.13) and X^ g by Eq.(3.15). 

3. By means of the measurement of the recoil electron spectrum n^^(T) in the ES 
process and the electron spectrum in the CC process. In fact, from Eqs.(3.29) and 
(4.2) we have 

$ = ^511 (45) 

X.,e(T)' ^^-'^ 

where E'^S(T) is given by Eq.(3.31) and X^^e(T) by Eq.(3.30). 

All these different methods of determination of the total flux $ will give results which 
must be in agreement with each other (otherwise the ratios Rnc) R^^(T), . . . considered in 
the previous section would be different from one). Let us stress that the proposed methods 
for the determination of the flux of ^B neutrinos do not depend on the mechanism of 
transition of solar z/g's into other active states. It is evident that a comparison of the 
flux $ determined directly from the experimental data with the ®B flux predicted by solar 
models will be an important test of these models. 

The SNO and S-K experiments will allow also to determine in a model independent 
way the probability of u^s to survive, Pi^^^,^^(E). In fact, we have 

P.^.(E) = 1^ , (4.6) 

where (p^^E) is the flux of solar z/e's on the earth and the total flux $ is given by Eqs.(4.3), 
(4.4) and (4.5). If it will occur that the probability Pj,^_>j,^(E) is less than one, this will be 
a model independent proof that z/g's transform into other states. A detailed investigation 
of the dependence on E of the probability Vy^^^^i^) will allow to distinguish different 
mechanisms of neutrino transitions (MSW, just-so vacuum oscillations and others) and to 
determine the corresponding parameters. 

The detection of solar neutrinos through the observation of ve -^ ve scattering will 
allow to obtain additional informations about the z/g survival probability. Let us define the 
average value of the probability of z/g to survive (Pi/e-+i/o)ES ^^ 

(P.e^.e)ES = y ^ y / K.e(^) " ^v,e{E)\ X(i?)P,,^,jE) dE , (4.7) 



^UeC ^^Vu.e 



-^th 



where 



X,,g= / a,^,{E)X{E)dE {U 



^T 
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^e -^ Vii{t) 



Am"^ (eV 



sin^ 26 



{^l^c^l^c) 



ES 



SMALL MIX MSW 



LARGE MIX MSW 



VACUUM OSC. 



6.1 X 10-6 
9.4 X 10-6 
8.0 X 10-" 



6.5 X 10-3 
0.62 
0.80 



0.32 
0.19 
0.31 



Table 4: Results of the calculations of the quantity {Puc~*uc)es ^^ ^^^ simplest model with i/g- 
i^nM mixing. The values of Am? and sin^ 29 used are also given. These values were obtained 
from the analysis of the existing experimental data (Ref.[20] for the MSW transitions and Ref.[21] 
for the vacuum oscillations). 

and Xj,^e is given by Eq.(3.15) (for E^ = 4.74 MeV, which corresponds to a kinetic energy 
threshold Tth = 4.5 MeV for the recoil electrons, we have X^^e = 2.12 x 10~^^cni^ and 
X^^e = 3.23 X 10"^^ cm^). 

The quantity (Pi/^^uJes ^^ connected with the total number of ES events A^^^ and the 
total flux $ by the following relation: 

-^ES 
Vvc-^i^c/ES 



Y Y 



$ 



-X. 



(4.9) 



The right-hand side of this relation contains only measurable and known quantities (the 
total flux $ is given by Eqs.(4.3), (4.4) and (4.5)). If the right-hand side of Eq.(4.9) will 
be found to have a value less than one, it will be a model independent proof that solar z/g's 
transform into other active states. 

Finally, a measurement of the recoil electron spectrum n^^(T) will allow to determine 
the average 



where 



1 



Xvee{T^) - X^^c(T) Je^{T) L 



^^-=^(i^,T)-^^(E,T) 



dT 



do-,. 



Em(T) 



dT 



dT 



-(E,T)X(E)dE 



X(E)P,^_.,jE)dE, 
(4.10) 

(4.11) 



and X,,^e(T) is given by Eq.(3.30). The functions X,,^e(T) and X,^^e(T) are plotted in Fig. 4. 
From Eq.(3.28) we easily obtain 



XUV - X.,e(T) 



$ 



XuAV 



(4.12) 
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Figure 10: Results of the calculations of tlie 
probability of fg's to survive, Puc^uc(^)i ^^ ^ 
function of the neutrino energy E, in the model 
with Ve'T^fiir) niixiug. The curves correspond to 
small and large mixing angle MSW transitions 
and to vacuum oscillations. The values of the 
parameters Am^ and sin 29 used in the calcu- 
lation are given in Table 4. 



Figure 11: Results of the calculation of 
(P uc^uc) -Es-T (^^^ Eq.(4.10)) in the model with 
Ve^i^fiir) mixing (T is the kinetic energy of the 
recoil electrons in the ES process). The curves 
correspond to small and large mixing angle 
MSW transitions and to vacuum oscillations. 
The values of the parameters Am? and sin^ 29 
used in the calculation are given in Table 4. 



If the right-hand side of Eq.(4.12) will be found to have a dependence on T, it will be a 
model independent proof that there are transitions of solar z/g's into other active neutrinos. 
We have calculated the survival probability Pj,^_»j,^(E) and the quantities (Py^-^u^)^^ 
and (P;yc^!^c)EST ^^ ^^^ model with Ve^^nir) mixing. The values of the parameters Am^ 
and sin^ 26 and the values of the quantity (Pj/c^!^c)es ^^^ given in Table 4. The values of the 
mixing parameters were obtained from a fit of the existing solar neutrino data under the 
assumption that the neutrino fluxes are given by the SSM. The results of the calculations 
of Pj,^_,,^^(E) and (Pi/o->i'o)est ^^^ presented in Figs. 10 and 11. These figures illustrate 
the possibility to distinguish the different mechanisms of neutrino transitions through the 
investigation of the energy dependence of the probability of z/g's to survive, Pj,^_^,y^(E), and 
of the averaged probability (Pj,^_,j,^)gg.rp. 



5 Conclusions 

We have shown that the existing solar neutrino data allow to exclude in a model indepen- 
dent way large regions of the values of the parameters Am^ and sin^ 29. We have also shown 
that from future solar neutrino experiments (SNO, Super-Kamiokande, Icarus) it will be 
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possible: 1) To reveal in a solar model independent way the presence of sterile neutrinos 
in the flux of solar neutrinos on the earth; 2) To determine directly from the experimental 
data the flux of initial ®B z/g's; 3) To obtain directly from the data the probability of Ve 
to survive. Thus, the future solar neutrino experiments will allow to check the predictions 
of the SSM independently from the properties of neutrinos and will allow to investigate 
neutrino mixing in a model independent way. 
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